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ABSTRACT 

Calculations have been made of fall-back disk heating by the pulsar wind as distinct 
from the soft X-rays emitted by the neutron-star surface. The relation between these 
heating rates and measured near-infrared fluxes in the K and Kg-bands places severe 
constraints on the inner radii of any fall-back disks that may be present in radio pulsars 
and in some anomalous X-ray pulsars. The lower limits found are so large that the 
disks concerned can have no significant effect on pulsar spin down. 
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1 INTRODUCTION 

Dust and debris disks are a commonplace phenomenon in 
astrophysics, being a feature of the major planets and of 
pre-main-sequence stars (see, for example, Beckwith et al 
1990). The properties of dust disks or of viscous ionized disks 
formed as a consequence of fall-back in supernovae have been 
a matter of continuing interest in neutron-star physics. But 
until recently, there has been no observational evidence that 
they exist. For a bibliography, we refer to Perna, Hernquist & 
Narayan (2000), Menou, Perna & Hernquist (2001), Black- 
man & Perna (2004), Ek§i & Alpar (2005), Ek§i, Hernquist 
& Narayan (2005) and to the many earlier papers on the 
subject cited by these authors. 

Further interest has been stimulated by flux measure- 
ments on the anomalous X-ray pulsar 4U 0142-1-61 in the 4.5 
and 8.0 /im infrared bands (Wang, Chakrabarty & Kaplan 

2006) . The energy spectrum derived from these measure- 
ments, and from earlier Ks-band and optical fluxes, is con- 
sistent with blackbody emission from a dust disk with an in- 
ner radius surface temperature of 1200 K, superimposed on 
a power-law spectrum. The optical emission, which is pulsed 
at the neutron-star rotation frequency, consists almost en- 
tirely of the power-law component, which is assumed to be 
of magnetospheric origin. 

There are also flux limits at 4.5, 8.0 and 24 yum in the 
mid-infrared for a limited number of isolated neutron stars 
(Wang, Kaspi & Higdon 2007) and for some supernova rem- 
nants at 4.5 and 8.0 pim (Wang, Kaplan & Chakrabarty 

2007) . But in the K and Ks-bands, a considerable number 
of magnitudes have been measured. These have been sum- 
marized by Mignani et al (2007) and include radio pulsars, 
anomalous X-ray pulsars (AXP), and soft-gamma repeaters 
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(SGR) . The purpose of the present paper is to see how these 
measurements constrain possible fall-back disk parameters. 

Fall-back disk formation is a detail of the supernova 
event and no attempt is made to describe it here. Initially, 
the disk must have been internally ionized and hence vis- 
cous, with mass and angular momentum transfer rates in- 
trinsic to its formation. We assume that at some later time, 
a transition to a passive disk occurs through the the ther- 
mal ionization instability (see Menou et al 2001) in which 
free electron recombination produces a sudden decrease in 
opacity which propagates rapidly inward. The disk formed 
has little internal ionization and viscosity, and is gaseous. 
Its luminosity is then principally a result of interaction with 
the pulsar wind or with blackbody radiation from the neu- 
tron star surface. The distribution of mass between dust 
grains and gaseous molecules depends on the various sub- 
limation temperatures concerned. Our assumption is that, 
at the time of observation, the inner radius of the disk lies 
outside the light-cylinder radius Rlc of the neutron star. 
An important consequence is that we can treat the disk as 
a partial termination of the pulsar wind rather than as a 
boundary condition to be satisfied by the solution for the 
magnetospheric fields. 

A previous paper (Jones 2007) on the interaction of a 
thin passive disk with the pulsar wind described the abla- 
tion processes in some detail and calculated the alignment 
and precession torques acting on the disk. Section 2 of the 
present paper extends the discussion of disk luminosity. The 
K and Ks-band magnitudes listed by Mignani et al for sets of 
radio pulsars and AXP or SGR are used, in Section 3, to con- 
strain the inner radii and temperatures of passive fall-back 
disks that might be present in these objects. Where pos- 
sible, the predicted mid-infrared fluxes from these objects 
are compared with the flux limits found by Wang, Kaspi & 
Higdon. 



2 P. B. Jones 



2 INTERACTION WITH THE PULSAR WIND 

The Deutsch vacuum solution for the electromagnetic fields 
E and B of a rotating neutron star (see Michel & Li 1999, 
Ek§i & Alpax 2005) makes it possible to write down explicit 
expressions for the spherical polar components of the mo- 
mentum density p, 
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in the inertial frame at radii r > Rlc, the z-axis being 
parallel with the neutron-star spin angular velocity fl. The 
momentum density in a physical pulsar wind, in which both 
fields and rclativistic particles are present, may differ from 
this in detail, but the time-averaged long-range part of the 
azimuthal component (p^) must be finite because it is di- 
rectly related to the neutron-star spin-down torque (see Sec- 
tion 2.1 of Jones 2007). The composition of the wind does 
not significantly affect its interaction with the disk and it is 
therefore correct to use this quantity derived from the vac- 
uum solution in calculations of ablation and, in particular, 
of alignment and precession torques. 

The vacuum-solution time-averaged polar component is 
(Pe) = for all angles, but its uni-directional time average, 
taken in the direction of either increasing or decreasing 9, 
can contribute to heating of the disk surface though not to a 
torque. With the addition of this component, the power in- 
put per unit area, to the one side of a thin disk that is viewed 
by an observer, can be written down and equated with the 
blackbody radiation emitted at the local temperature T{r), 
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where Cw is the pulsar wind luminosity, (3 is the disk tilt an- 
gle and ^ is the angle between the neutron-star magnetic and 
spin axes. The first term in equation (2) is derived from the 
long-range part of (p^) by means of equation (14) of Jones 
(2007). The second term, in the approximation of small /3, is 
from the long-range part of the uni-directional time average 
of pe evaluated in the vicinity old = 7r/2. The status of the 
vacuum solution form of this term is much less certain. Apart 
from an apparent unphysical singularity at ^ = 0, which is 
a consequence of the vacuum field solution and our use of 
Cw oc sin^ ^ as a working parameter, the term is linearly 
dependent on the radial component of the electric displace- 
ment D which is undetermined because the total charge of 
the star and the physical, as opposed to vacuum, magneto- 
sphere is an unknown quantity, as has been emphasized by 
Michel & Li (1999). 

The values that might be found for the tilt angle 13 
are, of course, unknown, but there is a case for assuming 
that its value at formation is non-zero as a consequence of 
the fact that type II supernovae are generally asymmetric 
(see Goldreich, Lai & Sahrling 1997). This suggests that the 
angular momenta of the neutron star, the fall-back disk and 
the ejected matter may not all be precisely parallel. 

The uncertainties in j3 and in the polar term are so 
considerable that we shall replace the vacuum-field quantity 
sin/3 -|- 2/3 cot ^ in equation (2) by a free parameter C, and 
adopt a conservative range of values 0.03 ^5 (" ^5 0-3- Ne- 
glect of the wind albedo is a further uncertainty in equation 
(2), but not one that is serious because relativistic parti- 
cles either in the wind or accelerated in the diffuse surface 



of the disk as described previously (Jones 2007) are atterm- 
ated only at depths of 10^ — 10^ g cm~^ within the disk. (We 
assume that the mean disk density E exceeds 10^ g cm~^.) 
For the same reason, it is not anticipated that infrared or 
optical emission from an undistorted disk would have an 
observable component pulsed at the neutron-star rotation 
period. 

We anticipate that, even for young neutron staxs, al- 
most all of the potential disk ablation will have occurred 

by the time of observation. Interaction with the pulsar wind 
increases the inner radius ri which is given by 
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which is valid provided the solution ri{t) exceeds the real 
Ti at formation. Here f2o is the neutron-star spin angular 
velocity at formation, I is the neutron-star moment of inertia 
and ruH the mass of the hydrogen atom. The remaining 
parameter is the number of baryons (mostly neutrons) 
removed from the disk per unit incident wind energy (see 
Jones 2007; equations 17 and 18). The mean mass density 
E is the initial value at formation. Owing to the number 
of unknown parameters it contains, equation (3) does not 
provide a useful constraint on disk inner radii derived from 
equation (2). But it shows that, in most cases, the existing 
ri must be near its f2 = limit. Because equation (3) is 
based on the interaction of the disk with the pulsar wind, it 
ceases to be valid when the light-cylinder radius exceeds r;. 
But for the radio pulsars, this occurs only for O so small as 
to be of no interest. 



3 DISK PARAMETERS 

Recent observations in the near-infrared have produced H, K 
and Ks-band magnitudes for many isolated neutron stars, in- 
cluding both radio pulsars and AXP. A useful listing of these 
has been published by Mignani et al (2007) who also in- 
clude estimates of the distances d. We have used the Spitzer 
Science Centre magnitude to flux-density converter to ob- 
tain the observed energy fluxes per unit frequency interval 
fu- Equation (2) provides the disk surface temperature for 
each object as a function of radius and it is then straight- 
forward, by integration and interpolation and with a num- 
ber of moderate assumptions, to flnd the disk inner radius 
that gives fv. In particular, we assume that the angle 
between the disk normal and the line of sight is fixed and 
given by cos^x = 0.5. No attempt has been made to cor- 
rect for the small K-band interstellar extinction. The values 
found for ri arc quite insensitive to the outer radius ro be- 
cause, for the temperatures found here, the outer part of 
the disk contributes little to the K-band flux. Therefore, we 
have adopted a fixed ro = 3.4ri, equal to the best-fit ratio 
found by Wang et al in the case of 4U 0142+61. With these 
assumptions, values of ri and the inner radius temperature 
Ti have been computed and are given in the right-hand six 
columns of Table 1 for C = 0.03, 0.1, 0.3. 

The disks of Table 1 need not be massive: for C, = 0.03, 
the masses are in the interval 10"'^ - lO"'"^ Mc5, assuming a 
density E = 10^ g cm~^, but their angular momenta can be 



of the same order as that of the neutron star. Certainly for 
the lower temperatures given in Table 1, the disks must con- 
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sist almost entirely of dust grains and only a small gaseous 
component. But even a very low concentration of free elec- 
trons and ions is sufHcient to stop the pulsar wind (see Jones 
2007) and there can be no doubt that such concentrations 
must bo present from a variety of sources. 

For consistency with our Section 1 definition of a fall- 
back disk, we require that a satisfactory solution has Vi > 
Rlc- The first six rows in the Table refer to radio pulsaxs, 
and they satisfy this condition by large margins. But it is 
interesting to sec that, in many cases, the AXP and SGR 
do not. In particular, there is no solution for 4U 0142+61, 
the object for which evidence of a disk is most compelling. 
This is not a cause for concern because equation (2) contains 
only the pulsar wind contribution to disk heating. For the 
AXP and SGR, the soft X-ray luminosity Cboi, possibly of 
thermal origin, exceeds jCw by several orders of magnitude 
so that its power input can be the more important if the 
disk is not of negligible thickness but subtends a finite solid 
angle at the neutron star, as in the model adopted by Wang 
et al (2006). 

It is reasonable to assume that a disk, immediately af- 
ter the ionization instability transition, would be plane and 
thin. If there is a gaseous component, it would have a diffuse 
surface of approximate depth 

h= IJ^, (4) 

where mA,z is an atomic (or molecular) mass and Ok is 
the Kepler angular frequency. For the 4U 0142-1-61 best-fit 
inner radius rt = 2.0 x 10^^ cm found by Wang ct al, the 
very small Kepler frequency, 1.5 x 10"'* rad s^^ leads to a 
fairly large depth, ft « 2 x 10^ cm, for atomic hydrogen at 
T = 10^ K. The 4U 0142-1-61 luminosity is jCboi ~ 10^^-10^®, 
so that although the disk is thin, the X-ray power input, 
(1 — rid)hCboi/ri exceeds that given by equation (2) even for 
the largo value of the X-ray albedo rjd adopted by Wang 
et al (2006). The solid angle subtended by the disk may 
also be increased by departures from the plane assumption 
arising from r-dcpcndcnt alignment torques (.loncs 2007) or 
disk instability (Pcttcrson 1977; Priuglc 1996). The X-ray 
contribution to disk heating is therefore important for the 
AXP and SGR, but not for the radio pulsars which have 

W ■ 



4 COMPARISON WITH OBSERVATION 

The electromagnetic spectra of a limited number of radio 

pulsars are known over wide frequency intervals. In the case 
of B0656-I-14, as an example, apart from the lack of measure- 
ments between W^^ and 10*'* Hz, the shape of the spectrum 
is of blackbody emission from the neutron-star surface su- 
perimposed on a background which decreases monotonically 
with frequency and is possibly of some common magneto- 
spheric origin (see Koptsevich et al 2001; Fig. 6) The K-band 
fiux listed in Table 1 is consistent with this spectrum. It is 
also known that the optical emission is pulsed at the rota- 
tion frequency. For example. Shearer et al (1997), observing 
B-band emission, found a la upper limit of only 8 x 10^'^* 
erg s~* cm~^ Hz~* on the unpulsed component of the flux. 
With regard to the AXP listed in Table 1, the optical emis- 
sion from 4U 0142-1-61 has a pulsed fraction of 0.3, and the 



optical and X-ray pulse profiles are in phase and of similar 
shape (Kern & Martin 2002, Dhillon et al 2005). 

Unambiguous evidence for the presence of a disk would 
be an unpulsed blackbody excess which, for the tempera- 
tures given in Table 1, would probably be best observed in 
the 4.5 and 8.0 /jm bands. With the possible exception of 
4U 0142-1-61, existing spectra provide no such evidence. 

Even for the smallest value of C assumed in Table 1, 
the radio pulsar disk inner radii that would be compatible 
with the K-band fluxes arc all larger than the best-fit ri — 
2.0 X 10" cm found by Wang et al (2006) in the case of 4U 
0142-1-61, although they are of a similar order of magnitude. 
It might be argued that the correct vahic of (" is much smaller 
than 0.03, perhaps by one or more orders of magnitude. For 
this to be so, the neutron star spin fl would have to be 
almost exactly normal to the plane of the disk. Also, the 
star and magnetosphere must be electrically neutral to the 
extent that the uni-directional time-average of pe, evaluated 
at 8 = 7r/2, is two or more orders of magnitude smaller 
than the vacuum solution value. Instead, we suggest that 
the inner radii given for the radio pulsars in Table 1 are 
broadly correct as lower limits. 

Heating of a disk with inner radius ri > Rlc by the 
pulsar wind is capable of producing the observed Ks-hand 
fluxes in only two of the AXP listed in Table 1. But these 
objects are differ significantly from radio pulsars in having 
very high soft X-ray luminosities, Cboi S> jCw, which give 
most of the disk heating, as in the model adopted by Wang et 
al (2006). The two exceptional AXP, J1048-5937 and J1810- 
197, are interesting in that 3a upper limits for their mid- 
infrared fluxes at 4.5, 8.0 and 24 ^m have been published 
recently by Wang, Kaspi & Higdon (2007) . These fluxes can 
also be calculated from the ri and Ti values given in the 
Table. For J1810-197, all the predicted fluxes are below the 
the 3(T upper limits. However, the ,11048-59.37 4.5 /um 3cr 
upper limit of 8 x 10^"^'' erg cm^'^ s^* Hz^* is exceeded 
by the predicted flux of 1.2 x 10^ erg cm" -'^ s"* Hz"* for 
( = 0.03 which increases to 3.1 x 10"^** erg cm"^ s"* Hz"* 
at C = 0.3. A very small value, C = 0.012 is required to 
fit both the measured Ks-band flux and the published 3a 
upper limit. In this case, the disk is at a higher temperature 
Ti — 1300 K and has a smaller radius ri = 7.9 x 10*" cm. 
But an alternative possibility is that the source of both the 
Kg-band and 4.5 /um fluxes is magnetospheric rather than 
a heated disk, which could be confirmed by a search for a 
pulsed component. 



5 CONCLUSIONS 

The present paper considers heating of a thin passive disk 
by non-radial components of the pulsar wind and computes 
inner radii and temperatures that are compatible with the 
observed K and Ks-band fluxes for a number of radio pulsars 
and AXP. In the case of the radio pulsars, the inner radii arc 
two or more orders of magnitude larger than the Rlc and 
the disks concerned, if they exist, can have no significant 
effect on pulsar spin-down. The conservative, and probably 
incorrect, assumption made is that the observed K and Ks 
fiuxes have no pulsed component and are of thermal origin. 
It is also assumed that there is no internal disk viscosity 
and that heating by the radial flux of soft X-rays from the 
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Table 1. Disk inner radii ri that are consistent with measured K and Kg-band magnitudes have been calculated, under the 
assumptions made in Section 3, for a number of radio pulsars and AXP or SGR. Inner radii temperatures Tj are also given. For 
J0142+61, J1810-197, IE 2259+586, and SGR 1806-20, the values of the wind luminosity Cyj have been obtained from the periods 
P and spin-down rates given, respectively, by Gavril & Kaspi (2002), Ibrahim et al (2004), Woods et al (2004) and Mereghetti et al 
(2005). For the remaining neutron stars, the periods and Cw arc from the Australia Telescope National Facility Pulsar Catalogue 
(Manchester ct al 2005). The distances d are those listed by Mignani et al (2007). The energy flux densities have been obtained 
from the observed K or Ks-band magnitudes using the Spitzer Science Centre magnitude to flux-density converter. The right-hand 
six columns give and Tj for f = 0.03, 0.1, 0.3. Blank spaces indicate the absence of a solution with r, > Rlc- 
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neutron-star surface is negligible except in the case of the 
AXP. 



The possibility that some of the individual neutron stars 
listed in Table 1 have disks has been investigated by previous 
authors. K-band fluxes from disks with inner radii near the 
Alfvcn surface have been calculated in the case of B0656+14, 
and the four AXP J0142+61, J1048-5937, J1708-4009 and 
IE 2259+586, by Perna, Hernquist & Narayan (2000). Also, 
the radio pulsars B0531+21 and B0833-45 were considered 
by Blackman & Perna (2004). The disk model differs from 
that of Section 2 in that the sources of heating were re- 
stricted to internal viscous dissipation and the soft X-ray 
spectrum of the neutron-star surface. Even so, the predicted 
K-band fluxes for the four AXP were ~ 10-2'^ - 10-2'^ erg 
cm-2 g-i Hz-i, considerably higher than the observed val- 
ues listed in Table 1. In the case of the three radio pulsars, 
inclusion of heating by non-radial wind components is cru- 
cial and leads us to the conclusion that the proposed disks 
are not consistent with the observed K-band fluxes. This is 
also true in the general case, treated by Ek§i & Alpar (2005), 
of disks with inner radii near the light-cylinder radius. This 
conclusion should perhaps be qualified because we do not 
consider disks with n < Rlc- Thus a disk satisfying this 
condition at formation, with Rlc given by the initial angu- 
lar frequency fio, would have its middle section removed by 
wind ablation leaving an outer section with the properties 
given in Table 1 and an inner section completely within the 
present light-cylinder radius. This inner sector would form 
part of the closed magnctosplicre and would not intersect 
open magnetic flux lines or interact with the pulsar wind. 
The properties of such remnant disks, if they exist, are not 
considered in this paper. 
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